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Abstract 

We obtain by dimensional reduction a (1 + 1) supersymmetric system introduced 
in the description of ultracold quantum gases. The correct supercharges are identified 
and their algebra is constructed. Finally we find static self-dual solutions that satisfy 
a Liouville type differential equation. 

Keywords: Super symmetry, Chern-Simons gauge theory, Ultracold gases, differen- 
tial equations, Liouville 

PACS numbers:11.30.Pb, 12.60.Jv, 67.856. Pq, 11.15.Yc, 02.30.Jr 



1 Introduction 



Since its birth in the early 1970's in the context of high energy physics and mathe- 
matical physics, supersymmetry has found a growing number of applications strongly 
influencing many areas both in experimental and theoretical physics. 

Originally proposed as a graded extension of the Poincare algebraJT], it was soon 
recognized that it can also be considered in the systems that exhibit Galilean invariance 
and this lead to the construction of a graded super-Galilean algebra[2] in d = 3 + 1 
space time dimensions. Constructing the supersymmetric extension of the Galilean 
invariant d = 2 + 1 Jackiw-Pi model[3], Leblanc et al [I] discovered the existence of 
2 graded superalgebras and related this to the possibility of finding BPS equations 
in the bosonic sector. After that, there were several developments related to Galilean 
supersymmetry in diverse contexts |5j. 

In (1 + 1) dimension Galilean supersymmetry was considered to study ultra cold 
quantum gasesjUElIH]. Ultracold quantum gases not only are interesting by the physics 
that they describe, but also are a useful tool in the modeling of others branches in 
physics [9j. One interesting example of this modeling was considered in Ref.[6, 7J. 
There, the authors propose that by combining a vortex line in a one-dimensional optical 
lattice with a fermionic gas bound to the vortex core, it is possible to tune the laser 
parameters such that a nonrelativistic supersymmetric string is created. This could 
allow to test experimentally several aspects of superstring and supersymmetry theory. 

From the theoretical point of view, the model that describe that proposal presents a 
supersymmetric structure. Despite the theory has interactions, the authors only found 
the generators of the super-Galilei algebra for the free theory. This fact constrast with 
a basic feature of supersymmetry theory, where the full Hamiltonian is generated by 
the supercharge algebra. Later, in ref.[8] was shown the existence of supersymmetry 
charge whose algebra generate the Hamiltonian with quartic interactions. Nevertheless 
a supersymmetry transformation associated to a charge be able to generate the full 
Hamiltonian of the theory was not found. In this paper we will show the existence 
of a supersymmetry related to a charge generating the full Hamiltonian of the theory 
presents on Ref.[6l [7]. Also we will construct the complete supersymmetry algebra. 
In order to discus this aspects we will present the (1 + 1) model study in Ref.[7] as 
the dimensional reduction of a Maxwell-Chern-Simons model proposed by Manton[10j. 
This has no influence in the derivation of the correct supersymmetry generators but 
leads us to the presences of interesting solitonic equations in the system. For the 
bosonic sector these equations reduce to a Liouville type differential equation. We 
analyze this case and construe the solucion for this type of Liouville equation. 

2 The Model 

Let us start by considering some features of the model proposed by Manton. This model 
is governed by a (2+l)-dimensional action consisting on a mixture from the standard 
Landau-Ginzburg and the Chern-Simons model, where the matter is represented by a 
complex scalar field (ft(x), 
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% + i) = j d 3 x( - X -B 2 + ij^d^ + ^ o |0| 2 ) - + 

k(A B + A 2 d A 1 ) + 1 A Q + A(|0| 2 - l) 2 - AiJf^j (1) 

Here 7, k and A are real constants, = da + iA^ (/i = 0,1,2) is the covariant 
derivative and B = d\A 2 — d±A 2 the magnetic field. The term 7^0 is related to 
the possibility of a condensate in the ground statepT] and jf is a constant transport 
current. It was show by Manton in Ref.[I.0J that this theory presents Galilean invariance 
with the requirement that the transport current transforms as jf — > Jj + jVi under 
a boost. With this consideration, we can choose a frame where jf = 0. The field 
equations in this frame takes the form 

i 7 A)0 = ~DiDi<f> - 2A(|0| 2 - 1)0 

€ijdjB = Ji + K€ijEj 

kB = 7 (H 2 - 1) (2) 
where Ei = OiAq — doAi is the electric field and Ji is the supercurrent defined by 

Ji = -I(^ A 0_0( A 0)t) (3) 

The first equation of this system is the non-linear Schrodinger equation. The second 
is the Ampere's law in two dimensions. The last equation is the Chern-Simons version 
of the Gauss law, which here takes a different form from that presents in the Jackiw-Pi 
model [3]. The energy of the system for static field configurations reads as 

E = J d 3 x(±B 2 + ±.( Dl ^D t( f> - \(\<P\ 2 - l) 2 ) (4) 

For finiteness we require that the energy vanishes asymptotically. This fixes the 
asymptotic behavior of the fields 

lim <f>(x) = e iaW , lim A t = d { a (5) 
where a is common phase angle. With these conditions the magnetic flux reads 

$ = f d 2 xB = f w=oo A i dx i = 2nN (6) 

where TV is a topological invariant which takes only integer values. Following to 
Hassaine et al.[12], we can rewrite the expression (|1J) as 
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E = / ± iD ^ 2 + + i - A) «i 2 - 1)2 T (7) 

where we have used the Gauss law of the equations ([2]) and the identity |-Di</>| 2 = 
\(Di±iD2)(f)\ 2z FB\(j)\ 2 ±e l idiJj. For A = =F 2^1 + 2"^ the potential terms cancel, and we 
see that the energy is bounded below by a multiple of the magnitude of the magnetic 
flux (for positive flux we choose the lower signs, and for negative flux we choose the 
upper signs): 

E>— 101 (8) 
~ 2m 1 1 v ; 

So this bound is saturated by fields obeying the self-duality equations 



(D 1 ±iD 2 )<P = 

kB = 7 M 2 -1) (9) 

Motivated by these results and the previous works on the (l+l)-dimensional supersym- 
metry in ultracold quantum gases [TJ El [8] we are interested here on the dimensional 
reduction of the supersymmetric extension of the model (pQ). Such extension can be 
carried out by considering the inclusion of non-relativistic (down-spinor) fermion V'[l] : 

5( 2+ i) = / d 3 x( - X -B 2 + ij^dtcj) + ^d t i, + iAM 2 + |Vf ]) - ^(A0) t A0 - 

^-{D^D^ + K (A B + A 2 d A l ) + 7 A) - ^-^Bi}j + Ai(H 2 - l) 2 + A 2 (|0| 2 - l)|Vf ) (10) 
2m 2m / 

where the coupling constants are given by 

2mn 2k z 2mn n z 

and we have include a Pauli term for the fermion corresponding to a down-spinor. 
This action is invariant under the following supersymmetry transformation 

5\(j) = V2mr]\ip , Slip = —\p2mx\\§ , 

(5iA = 0, <M° = - 7 = — (r^^t -ni4>^) • (12) 

if the coupling constants satisfy 

^-+2^-^=0 (13) 
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Where 771 appearing in (|12p is a complex Grassmann variable. In order to analyze 
the lineal problem [8| 1 13j. it is natural to consider a dimensional reduction of the action 
(jTUJ) by suppressing dependence on the second spacial coordinate, renaming A y as B. 
Then, the action (I10H becomes 



S(i+i) = f d 2 x( - \(d x Bf + i^dt^ + ^d t ip + iA p) - 1-(D X ^D X <P - ^-B 2 p - 

^-(D X ^D X ^ + K(A d x B + BdoAi) + 7 A + X 1 (p b - l) 2 + X 2 ( Pb - (14) 

Where we have introduced the matter densities, 

p b = \(/)\ 2 , Pf = \i[j\ 2 , p = p b + p f (15) 
The Gauss law constraint for this action is 



a x B = l(fl-i) (i6) 

K 

Note that this constraint has an additional constant term from that appearing on 
Ref.([Hl [13]). The equation can be solved as 



B{x) = i J dz < x ~ z )(p{ z )- 1 ) ( l7 ) 

Using these expressions for the magnetic field and its derivative, the action ([1 
takes the form 



S {1+1) = f d 2 x(i 7 (^d t cP + V^V) - ^(D x <j>?D x <f> - ^(D x ^D x tP - ^-B 2 p + 

Here the action contains new coupling constants defined as 

2 2 

A 'i = A i - h ' A 2 = A2 - 'h ~ v- ( 19 ) 

Following the method proposed in Ref.|13|. the gauge field A x may be eliminate 
from the action (|18p via a gauge transformation. Indeed, after transforming the matter 
fields as 

<j){x) ^ e- ia{x) 4>{x), il)(x) -»• e- ia ^iP(x) (20) 

with 
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2 

the action can be written simply as 



a(x) = - / dze{x - z)A x (z) (21) 



\' 1 (p b -l) 2 + \' 2 (p b -l)p f -l-B 2 p) (22) 
The last term appearing in this action is a constant of motion, 

1 ( ^ f 

— / dxiB 2 (xi)p(xi) = g 2 \ dx\dx 2 dx 3 e(x\ - x 2 )e(x\ - x 3 )(p(x\) - \) 



7 2 f 

(p{x 2 ) - l)p(x 3 ) = 24 mK 2 / dx 1 dx 2 dx 3 (p(x 1 ) - l){p(x 2 ) - l)p(x 3 ) 

(N 3 - 2N 2 + iV) 7 2 
24mK 2 

where N = J dxp(x), and use has been made of the identity 



(23) 



e(xi - x 2 )e(x 1 - x 3 ) + c(x 2 - x 3 )e(x 2 - x%) + e(x 3 - x 1 )e(x 3 - x 2 ) = 1 . (24) 
Thus, dropping this term, the action can be written as 



i+i) = J d 2 x ( i7(0t$0 + ^d t ij) - l-{d x ^d x 4> - ^-(d x 0d x ^ + 

A^-l^ + A^-l)?/) (25) 



This is the model study in Ref. [6] , in the context of description of ultracold gases 
dynamics. The action (|25p is similar to the action derived by dimensional reduction 
of the Jackiw-Pi model [8]. The difference is that our model also include the chemical 
potential terms. 

Let us now consider the possible supersymmetry transformations that leave un- 
change the action (|25p . One obvious supersymmetry of the system (|25p takes place 
when bosons and fermions are interchanged according to 

5i4> = \/2mri\ip , Slip = — \/2mrji4> , (26) 
and the coupling constants are related by 

\ 2 = 2A' X (27) 



6 



It is interesting to note that the transformation (|26p is also a super symmetry of the 
model study in Ref . |8j . 

The fact that is less evident is the existence of a second super symmetry. In dis- 
cussing this super symmetry, note that the the action (f25l) is invariant under the fol- 
lowing transformation 

h<f> = -^n\{dx^ - Bii)) + ib 2 a (j> , S 2 if) = --^V2(d x 4> + Bi4>) + i5 2 a if) , 

(28) 

provided that 



a; 



7 
2mn 



A 2 — 2A X 



(29) 



and 



Bi{x) = / dze(x - z)p(z) 
2k 



5oa 
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2K\/2m 



dze(x — z)(r]24>ip^ — r/l 



(30) 



Notice that combining equations (|19p and (|29p we obtain Eq. ljlip . As the transfor- 
mation ()26p . the expression (I28p is also a supersymmetry of model explored in Ref.[8j. 
Another interesting fact is that the condition for coupling constants imposed in (|29p is 
a particular case of the equation (I27p . This implies that the following combination of 
the precedent two supersymmetries 



'2m 



rjlidxtp - Blip) + ib 2 a 4> - '-^j=rip 



5 2 if> 



--r/ 2 (d x (f) + Bi4>) + i5 2 a if) + 



/Km 



V2m ' V Km 
is also a supersymmetry of the accion (|25p if the condition (|29p is kept. 



(31) 



3 The supersymmetry algebra 

In this section we shall study the algebra of the generators associate to the transfor- 
mations ([26D and 

In the representation of the supersymmetry algebra the generator associated to the 
supersymmetry (|26|) can easily be fond to be 

Qi = -iV2m J dxip ] (p (32) 

In order to write the supersymmetry algebra we define the Poisson brackets for the 
functions of the matter field as 



7 



r pr n -■ [ a ( 6F 5G 5F SG 5TP 51 G Sr ' F 51 G \ (v* 

1 ' ]PB l J \5ft(r)5<j>(r) 5<f>{r) 8ft (r) S^(r)S^(r) #(r) 5^{r) ) 1 ' 

where the subscripts r and I refer to right and left derivatives and in particular we 
have 



{4>(x 1 ,t),(f)*(x2,t)} = -i5(xi - x 2 ) {ip(x 1 ,t),ip*(x 2 ,t)} = -i5(xi - x 2 ) (34) 
Using the definition of Poisson bracket it is easy to get 

{Q 1; Q\} = -2im j dxp = -2iM . (35) 
Following the Ref.[7j we can define a second generator 

R = tL= f dx^d x 4> . (36) 

y2m J 

which generate the free part of the Hamiltonian 

{R, R^} = ^— J dx {ftd 2 x <P + ^d 2 x ^) = -iH free (37) 

Nevertheless, as we mentioned in the introduction, the previous works does not 
propose a supercharge be able to generate the full Hamiltonian of the model explored 
in Ref.[7], that is the Hamiltonian derived from the action (I25p 



H 



d 2 x(^-{d x <trfd x <t> + -^-{d x ^d x i> - 
\2m 2m 

1 iPb-if-^ziPb-VPf) (38) 



2m/t rriK 

In addition we know, from the Ref.[8], that the supersymmetry (|28[) is generated 
by the supercharge 

Q 2 ] = i= [ dx^(d x + B 1 ) 



/2m 

= / dx^(x)(d x + ^- [ dze(x - z)p(z))<f>(x) (39) 

which generate a Hamiltonian only with potential terms quartic in the fields. Based 
in this idea and in the fact that the transformation (|3ip is a combination of the super- 
symmetry (f2"6"|) and (|2S|) . it seems natural to define a supercharge Q 2 , such that it is a 
linear combination of Q\ and Q2 



<s 



Q 2 = i= I dx^{d x + B x )<t>+^P= [dxifi<f> (40) 

V2m J \/Km J 

We will show that this charge generate the Hamiltonian (j38[) . Using this charge 
we can calculate 

{Q 2 ,QI} = {(Kx + K 2 ), (K\ + K$)} = 

{K U K{} + {^1,4} + {k 2 ,K\} + {k 2 ,K\} (41) 



where 



K X = Q^ K 2 = ^fLQ 1 (42) 



The first of the brackets was calculated in Ref. [8] and its result is 

k x ,kI) = ~ [d 2 x( J- (d x ^d x <(> + (d x ^d x ij 

) 2m J \2m 2m 



+ 

\'Zm ' ' 'lm ' 

B lP ~ -PfPb ~ d x B x p b + d x Bip f ) (43) 

which can be reduced, after using the definition of the B\ and eliminating the 
constant term J d 2 xB 2 p, to 

\k u kI} = [ d 2 x(^-(d x ^d x <p + ^-(d x ^d x ^- 

I J 2m 7 \2m 2m 

2 -P/^ " ( 4 4) 



K ft 

The second bracket can be developed as follows 



[k 1} KI} = ~^J^j dx l dx 2 ( K {^\x l )^{x 2 )}4>\x 2 ){d x + B 1 {x 1 ))4 ) {x l ) + 

tf(x x ) [d x <t>{x{), 4t(x 2 )} tf>(x 2 ) + ipHxJ {b 1 (x 1 ), tffaMxt)} 4>{xi) + 
^(x 1 )B 1 (x 1 ) {<f>( Xl ), <^(x 2 )} rP(x 2 )) (45) 
It can be easily checked that 

{b 1 (x 1 ),^(x 2 )^(x 2 )} =0 (46) 

Thus 



\K u Kl} = [dxUi(x)(d s + B 1 (x))<l>(x) + 

J m\J2n J \ 

^{x){d x + B x {x))^{x)) 
In similar form we have for the third bracket 

\k 2 ,k\\ = l — [ dxU ] {x){d x - Bi(x))<t>(x) + 

•> ) m\J2n J \ 

From (@7J) and (gSJ) we get 

{tfi,tfj) + {tf 2j tf}) = —^=[dxB l (x)p = 

/dxdzeix — z)p(x)p(z) - 

The last bracket gives 



Then the full bracket takes the form 



27 7 2 , 2 7 \ ■ rr 

— P/Pfc - -ft + — p = -»•« 
The algebra is completed by the following bracket 



{qi,q\} = {qi,kI} + {q 1 ,kI} 

The first of this brackets can be calculated to give 



Where we have used that 



J dxBi(x)p = 



10 



The second bracket is 



dxp (55) 



Therefore 



i J d 2 x($d x (j) - d x $$ + ^d x ip - d 2 

1 dxp (56) 



In order to show that this expression may be identified with the linear momentum 
of the system we use the Noether's theorem. Let {9 C } = {<ft, (fv , tp, ip'} the set of the 
fields of the our system, where c runs from 1 to 4. The theorem establish that if under 
a variation of the fields 59 c , the variation of the Lagrangian density is a surface term, 
5C = d^X^, then exist a conserved current associated with such variation of the fields. 
The Noether current, assuming the summation convention over the index c, is 

J^wm se <- X " (57) 

We are interested on the zero component of this current associated to the transfor- 
mations 

{se c } = {d x <i>, o x <p\d x ip, d x ^} (58) 

where 



X° = -iyJ^-p, X 1 = 0, (59) 

Note that, from the Noether theorem, the only restriction for X^ is that d^X^ be 
a surface term. Using (|57p . the linear momentum is 



P = j fdx = -i 1 {Q 1 ,Q\] (60) 
The energy-momentum-tensor, then may be defined as 

T « = iJInr^ ~ C6 » + *7\/%^ (6i) 

0(0^0) V K 

where 6„ = 1 if p. ^ v and 5u = if p = v. 

Finally it is easy to check that the remaining brackets are zero 

{Qa,Qp} = {QlQl} =0 (62) 
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4 The soliton solution 



Consider now the derivation of the self-dual equations. As discussed in Ref. [14] the 
field B plays an important role in the derivation of self-dual equations. Indeed the 
expression (1171) of B involves the existence of a novel soliton. The action (122|) may be 
easily reexpressed as 

s = t SxU^d^ + ^dt^-^^ + cB)^ 2 -^\(d x + cB)^\ 2 - 

J \ 2m 2m 

^-d x Bp + >: i (p b - If + \' 2 ( Pb - l) Pf ) , (63) 



where £ = ±1. Using the Gauss law f|16f) and after a bit of algebra we have 

1 . .., 1 



(a; - - 1) 2 + (A' 2 - ^)( Pb - i )Pt - S^b) , (64) 

which leads, in the static field configuration, to the Hamiltonian of the form, 



H = Jd 2 x(^\(d x + tB)<t>\ 2 + ^-\(d x + (B)i>\ 2 - 

(A'i - - I) 2 - (A; - %(p 6 - l) P f + £-d x B) (65) 

2mn mn 2m / 



We can choose Ai = S- 1 — and An = — so that our Hamiltonian becomes 

l Iran £ ran 

H = l d2x (^ {dx+cm2+ ^ {dx+cB) ^ 2+ S^ dxB ) (66) 

The last integral vanish since B must be zero in de boundary. Then, at the minimum 
of the energy configurations, the self-dual equations are satisfied 



(8 X + (B)^ = (67) 
(d x + (B)i; = 



Notice that for the particular choice £ = 1, we recover the supersymmetric case. 
This equations can be explicitly written by using the equation (TjFi 



d x 4>{x) + ^ ( y dze ( x ~ z)p(z)(j)(x) - J dze(x - z)4>{xfj = (69) 
d x ip + p- ( ! dze(x - z)p(z)ip(x) - I dze{x - z)tp(x)j = (70) 
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which present an additional linear term from those found in Ref . [8] . When ip is set 
to zero, the above set of equations reduces to 



d x (j>{x) + ^ J dze(x - z)( Pb (z) - \)4>{x) = 
Assuming a solution of the form </> = y/pH, we arrive to 



(71) 



ld x (logp b (x)) + p- [ dze(x - z)( Pb (z) - 1) = (72) 
2 Ik J 

Differentiating with respect to x, we get the following one-dimensional Liouville 
type equation 



i^(logp 6 (x)) + ^(p 6 (x)-l) = 
We propose as the solution of this equation the following series 



(73) 



1 + OnSech n (bx) 



(74) 



n=l 



Here a n are the real coefficients of series, b is a real constant and we have renamed 
p as Pb{x). In order to check that this is really a solution we rewrite Eq. (f73l ) as follows 



(d xP y + (d 2 x p)p + vp 2 (p-i) = o 



where v = When the series (|74p is introduced into Eq. (|75p we obtain 



(75) 



oo 

s [~ anc 



a m mnb 2 + n 2 a n a m b z + 2va n a Tl 



n,m=l 
oo 



sech n+m (bx) + 
sech n+m+2 (bx) - 
^(n 2 + n)a n 6 2 sech 2+n (bx) + ^ [n 2 a n b 2 + v& Q sech n (bx) + 



\a n a m nmb 2 — n 2 a n a m b 2 — na n a m b 2 



n,m=l 

oo 



n=l 



n=l 



OO 



v anamOjsech n+m+i (bx) = 

n,i,m=l 



where for arriving to this expression we have used the relation 



(76) 



tanh 2 (6x) = 1 — sech 2 (bx) 



(77) 
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So we have an expansion of powers of sech(bx) which must be equal to zero. This 
implies that the coefficient of each power must vanish separately. Prom the coefficient 
of sech(bx) we obtain 



v 



(78) 



whereas from the coefficients of sech 2 (bx) and sech 3 



(bx) we have 




(79) 



and 



a-3 = 



3a? + 2ai 
8 



(80) 



The method can be continued in order to determine the rest of the coefficients. 

5 Conclusion 

In this article we have studied a (l+l)-dimensional model introduced in the descrip- 
tion of the supersymmetric-ultracold gases. This model and its super symmetries was 
previously studied in Ref . [6j [71 [8] . However the problem of finding the supersymmetry 
algebra that generate the full theory was unresolved. In this note we started by ex- 
tending supersymmetrically a model proposed by Manton, and related it to the theory 
of ultracold gases. Then, the correct supercharge that generate the full theory were 
identified and their algebra was constructed. In addition the solitonic structure was 
analyzed and novel solitons are found. 
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